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In this talk, we report our recent studies on an effective thermodynamic potential (f^cff) at fi- 
nite temperature (T 7^ 0) and zero quark-chemical potential (/ir — 0), using the singular-gauge 
instanton solution and Matsubara formula for Nc — 3 and Nf = 2 in the chiral hmit, i.e. rUq — 0. 
The momentum-dependent constituent-quark mass is computed as a function of T, together with 
the Harrington-Shepard caloron solution in the large-A'^c hmit. In addition, we take into account 
the imaginary quark-chemical potential = A4, indentified as the traced Polayakov-loop ($) as 
an order parameter for the Z{Nc) symmetry, characterizing the confinement (intact) and deconfine- 
ment (spontaneously broken) phases. As a consequence, we observe the crossover of the chiral (x) 
order parameter and It also turns out that the critical temperature for the deconfinement 
phase transition, is lowered by about (5 ~ 10)% in comparison to the case with the constant 
constituent-quark mass. This behavior can be understood by considerable effects from the partial 
chiral restoration and nontrivial QCD vacuum on the <1>. Numerical results show that the crossover 
transitions occur at {T^,T^) ~ (216,227) MeV. 
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I. INTRODUCTION 



We note that the phase structure of quantum chromodynamics (QCD), as a function of temperature T and quark- 
chemical potential fi, represents the breaking patterns of the relevant symmetries in QCD. Simultaneously, each 
QCD phase can be characterized by the corresponding order parameters, reflecting the nature of the symmetries. In 
this sense, exploring the QCD phase diagram is of great importance in understanding strongly interacting systems. 
Especially, recent energetic progresses, achieved in the ultra-high energy experimental facilities, such as the RHIC, 
have triggered much interest to investigate the QCD phase structure in the vicinity of high T « T^, whereas ^ remains 
relatively small, being similar to the early universe. 

Starting from the first principle, the lattice QCD (LQCD) simulations must be a promising method to investigate 
this region (T 7^ and ^ ~ 0) with less difficulties, such as the sign problem [IHS]. Many attempts have been also 
done in various effective field-theoretical approaches [THTO]. Among them, interestingly enough, the Polyakov-loop- 
augmented Nambu-Jona-Lasinio (pNJL) model describes the crossover of the two different QCD order parameters 
for the chiral and Z(iVc) symmetries, represented by the chiral condensate (qq) cx cr^ and the traced Polyakov loop 
{(j)) = respectively [T HTTl [H] . If the 'Z{Nc) symmetry is intact, $ becomes zero, indicating the confinement phase. 
On the contrary, provided that the symmetry is broken spontaneously, one has $ 7^ for the deconfinement one. 

Instanton model can be also thought as an appropriate framework to be employed for this finite-T subject, con- 
sidering that it has provided remarkable descriptions so far for various nonperturbative QCD and hadron properties. 
Note that the instanton solution at finite T, being periodic in Euclidean time, i.e. caloron turned out to be essential 
for this purpose [501 [H]- Nonetheless for its relevance, its practical application is still under development P^27j . 
Confinement properties have been discussed as well with semi-classical objects, such as the meron (a half of regular- 
gauge instanton), by indicating the area law for the Wilson loop [281 [29]. The caloron with non-trivial holonomy, so 
called the Kraan-van Baal-Lee-Lu (KvBLL) caloron [^HIISI], was taken into account as a lump of dyons to understand 
the QCD confinement 22, 23, 25 . 

In the present talk, we want to develope an effective thermodynamic potential (fioff) at finite T with /iR = 0, 
employing the instanton framework. Our strategy is rather simple and practical as follows: 

i) Using the instanton distribution function at finite T from the caloron solution with trivial holonomy (the 
Harrington-Shepard caloron) [301l31j . we first compute the instanton density and average size of instanton as functions 
of T, resulting in that the instanton effect remains finite even beyond the critical temperature Tc ~ Aqcd- Taking 
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into account these ingredients, we finally obtain (fe(tliree momentum), T(temperature))-dependent constituent-quark 
mass M, M^^t, which plays the most important role in the present approach. 

ii) In constructing r^effi "we take into account a practical way, instead of using the caloron and its quark zero-mode 
solution: ficfi is obtained by applying the Matsubara formula to the effective action, which is derived from the usual 
singular-gauge instanton solution at T = 0, as done usually in effective models [7 HTT1 [TlHl9j . 

iii) The singular-gauge instanton solution is nothing to do with the confinement |29j . On the contrary, it explains the 
nonperturbative QCD properties very well in terms of the spontaneous breakdown of chiral symmetry (SB^S). Hence, 
considering the chiral and Z{Nc) symmetries on the same footing as in the pNJL model, we introduce the imaginary 
quark-chemical potential ni = A4, which corresponds to the uniform color gauge field in the Polyakov gauge. It will 
be indentified later as the traced Polyakov loop as an order parameter for the spontaneous breakdown of Z(iVc) 
symmetry, i.e. the deconfinement phase transition. 

Considering all these ingredients, now we can write a neat expression for QcS as a function of T with Mk^r- By 
solving the saddle-point equations with respect to the mean fields, i.e. a and we can draw the curves for ct^ and 
$ as functions of T. From the numerical calculations, we observe that the first-order deconfinment phase transition 
in pure-glue QCD is modified to the crossover one, according to the mixing of dynamical quarks and in ficff- In 
contrast, the mixing gives only negligible modifications on cr^. As a result, the crossover of the two different QCD order 
parameters is shown by the mixing. It also turns out that Tc for the deconfinment phase transition (T^) is lowered by 
about 10% with Mk.Ti in comparison to the case with a constant M (Mq.o ~ 350 MeV). From this observation, we can 
conclude that the nontrivial QCD vacuum contribution and partial chiral restoration play a considerable role even for 
deconfinment phase transition. The numerical results show that = 216 MeV and Tf = 227 MeV. The discrepancy 
between them becomes about 10 MeV, which is rather larger than that computed in the local pNJL model [SUlOj. 
We note that the LQCD simulations provides smaller values than ours, whereas Tf is comparable [T]-[5] . We also find 
that depends much on the partial chiral restoration. 

We organize the present report as follows: In Section 2, we briefiy explain the theoretical framework for obtaining 
the effective thermodynamic potential. The numerical results with discussions are given in Section 3. The final 
Section is devoted to the summary and outlook. We note that this talk is based on Ref. , and more details on the 
theoretical evaluations can be found there. 



II. THEORETICAL FRAMEWORK 



In the present section, we drive an effective thermodynamic potential fioff , considering all the ingredients discussed 
in Section II and III of Ref. |32j. In addition, we take into account the imaginary quark chemical potential (/^i = ^4), 
which will be indentified as the traced Polyakov loop as an order parameter for the deconfinment phase transition [71 
[T4] . Moreover, this corresponds to the uniform color gauge field, induced in the Polyakov gauge, = (0,^4). All 
calculations will be performed in the case for = 3, Nj = 2, and /iR = in the chiral limit for the leading large- A^^, 
contributions. In order to evaluate ficff as a function of T from the effective action, we employ the anti-periodic 
Matsubara formula for fermions in Euclidean space: 
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where the Matsubara frequency reads w„ = (2n+l)7rr. We also include an imaginary quark-chemical potential [TlllOj. 
which can be identified as the fourth component of SU(A'c) gauge field {A4) in Euclidean space, resulting in a simple 
replacement /c — > fc — A4 in ficff- Using these ingredients, we can have the following effective thermodynamic potential 
per unit volume in the presence of the dynamical quarks and A4 : 
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where the cr is a function of T as well. Note that we have ignored A4 in M for simplicity. After a tedious but 
straightforward manipulation, one is led to the following expression: 
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FIG. 1: Mfe T in Eq. |5} as a function of T and the absolute value of three momentum |fc| [GeV]. 



Here, we have used that 
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where the Mq^t is the T-dependent constituent-quark mass for fc = 0. For more details, see also Appendix in Ref. |32j . 
The numerical result for M^ t is shown in Figure [l] as a function of T and |fe|. From the figure, it turns out that the 
fc dependence of M becomes weak as T increases. At T = 0.4 GeV, the fc dependence does not appear at all. This 
behavior can be understood as follows: The k dependence of M is generated from the quark-instanton interaction, 
i.e, the delocalization of the quark fields in the instanton ensemble. As the instanton ensemble becomes more dilute 
as T increases (smaller p, equivalently) , the interaction probability decreases, resulting in the reduction of the k 
dependence. However, note that the dependence still remains visible around Tc ^ Aqcd ~ 200 MeV as shown in the 
figure. 

Now, we are in a position to consider the traced Polyakov loop 0, defined in a SU(-/Vc) gauge group as: 



-Tr.exp — 



-Tr.exp — 



(6) 



Taking into account the Polyakov gauge, A4 is diagonal in a A'^ x Nc matrix. For instance, the perturbative YM 
potential can be expressed by this physical quantity [35] and prefers the deconfinment phase [101111]: (0) becomes 
finite according to the spontaneous breakdown of the Z^^ symmetry with the trivial holonomy. In contrast, if the 
symmetry is intact, one finds ((/)) = 0, indicating the confining phase with the non-trivial holonomy. 

Thus, plays the role of an exact order parameter for the 'L{Nc) symmetry for pure-glue QCD, in which the quark 
degree of freedom is decoupled according to its infinitely heavy mass. It is worth mentioning that the 'Z(Nc) symmetry 
is broken explicitly in the presence of dynamical quarks with finite mass, considering the anti-symmetric nature of 
fermions, resulting in that (0) is not an exact order parameter for Zj^^ symmetry any more. However, from the 
phenomenological point of view, incorporating (j) and dynamical quarks has been quite successful to a certain extent 
to explain various features of the QCD phases transition: The crossover near Tc for instance for Nf — 2. Hence, in 
the present talk as done in the pNJL model, we have incorporated the instanton-based model with the SB^S and the 
traced Polyakov loop as an order parameter for the Z(Ac) symmetry. The trace over color space in Eq. ^ can be 
evaluated further in terms of (/) and (j)* using Eq. ^ as follows: 
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On top of rJ^Jj^^" (now becoming fi^g"^), an additional pure-glue effective thermodynamic potential was suggested 
in Refs. [gHTTl [ftj as a function of (j) and 0*: 
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where the coefRcient 62 is a function of T: 
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Here, the Tq denotes critical T in pure-glue QCD, resulting in Tq = 270 MeV at which the first-order phase transition 
occurs. Note that Tq is different from Tc, which will be computed later in the presence of the mixing of the dynamical 
quarks and (j). The coefhcients, a and 6, are listed in Table |l] This parameterization of the effective potential 
in Eq. ([s]) bears the Z^^ symmetry, conserved in pure-glue QCD by construction, and works qualitatively well up to 
T sa (2 — 3) Tj,, from which the transverse gluons come into play considerably. 

Finally, substituting Eq. ([t]) into Eq. ([3| and adding Eq. ([s]) to it, we arrive at the following expression for the effective 
thermodynamic potential with the two order parameters, cr^ for the chiral phase and (j) {<p*) for the deconfinment 
phase transitions, at finite T and fiR = 0: 
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where we have replaced the (p into its mean value {(j)) = <1> as done in Refs. [StilOj. Although, this expression for fioff 
is very similar to those given in Refs. [8l41H [T4]. ours is distinctive from them quantitatively in several points: 

i) The scale parameter of the model Aw 1/p is obtained as a function of T by solving the instanton distribution 
function as discussed in the previous Section. Moreover, the M is expressed as a function of k and T (Mfe.T), rather 
than a constant, manifesting the partial chiral restoration and nontrivial QCD vacuum contrib utio ns. 

ii) Consequently, there appears no divergence in the energy integral oc J cfik E^^t in Eq. (10 1 by virtue of the 
fe-dcpendent M, which plays the role of an intrinsic ultraviolet (UV) regulator, being different from other local- 
interaction models, such as the usual pNJL model. 

iii) At the smae time, as for the 2Nf-t Hooft interaction, the quark-meson coupling strength depends on k as well 
as T, being different from that in other models, in which it is a constant value fixed at zero T. 

iv) All the relevant quantities at zero T, the Mo^o: and rig for instance, are determined self-consistently by solving 
the saddle-point equations. 

Now, we evaluate the equations of motion with respect to the mean fields, ct, $, and $, by minimizing fJeg: 
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From the first two equations, it can be easily seen that $ = at the saddle point, we can compute the values for the 
and $, resulting in the following asymptotic behaviors: 
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showing appropriate asymptotic chiral behaviors as expected. 
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TABLE I: Coefficients for 0.'^^^ in Eqs. |8f and |9|, taken from Refs. [MlOl 
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FIG. 2: "1> (increasing curves) and normalizaed cr^ (decreasing ones) for To = 270 MeV using AI^^t- The solid and dashed hnes 
indicate the cases with and without the dynamical quark and $ mixing, respectively. 



III. NUMERICAL RESULTS 



In this Section, we present the numerical results for the two different order parameters, cr^ and $ as functions of T, 
and related discussions are given as well. Hereafter, we will used the normalized value for cr^, ct^/cto for convenience. 
Firstly, in Figure [2j we draw them for the cases with (solid lines) and without (dashed lines) the mixing of the 
dynamical quark and $, using M^^t- Here, we employ Tg — 270 MeV for the pure-glue potential in Eq. ([S]). In order 
for the unmixed case (dashed lines), in which the quarks are decoupled from the pure giuodynamics, we set M infinite, 
equivalently Ek^T — > oo. As shown in the figure, even with or without the mixing, showed the crossover for the 
chiral phase transition. On the contrary, the phase transition pattern for ^ turns out to be distinctive, depending on 
the mixing. In this sense, the crossover of <I> is caused by the mixing, as suggested by the pNJL models [7HT01 [Tl] and 
shown by the LQCD simulations with dynamical quarks [B]. An interesting behavior shown in Figure [2] is that the 
quark-gluon mixing effect increases then decreases for tr^ as indicated in the solid and dashed lines. This tendency 
comes from the combination of the opposite behaviors of $ (increasing) and e"^*"-^/^ (decreasing) in a^. 

In Figure [3j we show the numerical results for $ as a function of T for different types of M, as listed in Table [III We 
again employed Tq — 270 MeV. $ for the pure-glue potential, showing the first-order deconfinment phase transition, 
is also given in Figure [s] for comparison. From the figure, it is clearly shown that the T- and/or fc-dependent M make 
<& shifted to lower T, resulting in lowering Tc (we will discuss the numerical values for Tc later in detail). 

For simplicity, we approximate it considering only the leading contributions as follows: 
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As known from Figure [T] and Eq. ([5|, when we take into account the T- and/or fe-dependent M, the strength of M 
decreases as k and/or T increases, and the same for Ek^x- As a result, the integrand in the r.h.s. of Eq. (13) tends 
to be larger, as k and/or T increases, than that with a constant M, Mo,o 
T = Tq in Ec^. ( 13 1 for example, then have 
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TABLE II; Notations for M given in Figure |3] 
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FIG. 3: <E> as a function of T for To = 270 MeV, employing Mk,T (solid), M^fi (dotted), Mq.t (dashed), and Mo,o (long-dashed). 
We also draw <t> for the pure-glue case (dot-dashed). 



where we have used the notation, Eqq = Mqq + k'^. In order to satisfy the relations of Eq. (14), $ for A, must 
be bigger than $b at T ~ Tq, since $ is positive definite and 64 ^ 63. This observation is also true for arbitrary T, 
except for the limiting cases, T — )■ or 00. Ats the same time, therefore, this situation can be understood as that $ 
is shifted downward almost horizontally: is lowered consequently. 

However, this lowering Tc behavior must be understood separately for the k and T dependences in M, since the 
both curves with Mo,t and Mk^ show it, as depicted in Figure [Sj 

i) T dependence: Mq^t 

It tells us that the partial chiral restoration, indicated by decreasing M ex fP with respect to T, has effects on the 
deconfinment phase transition, although these two phase-transition mechanisms are believed to be different to each 
other. Here is a microscopic explanation for this: As T increases, the pseudo-particle (instanton) become smaller in 
its size ~ yO, decreasing QCD vacuum contribution simultaneously, resulting in that the quarks are enable to travel 
more freely with less interactions with the instantons, and strings attached to each quark are extended more at lower 
string tension, in comparison to the case without the partial chiral restoration in M. Hence, the condensation of the 
strings can happen easily at lower T, toward the deconfinment phase. 

ii) k dependence: M^^o 

The k dependence in M is originated from the delocalization of quark fields in the presence of the instanton 
background [33^, not from T-related effects. As momentum transfer increases, quarks become lighter cx l/k*, loosing 
the nontrivial QCD vacuum contributions. In other words, the instanton effect is weakened seemingly in the k 
integrals. Consequently, being similar to the T-dependence case, the deconfinment phase transition takes place at 
lower T. 

From these explanations, one is led to a conclusion that the QCD vacuum contributions and the partial chiral 
restoration play considerable roles for the deconfinment phase transition to a certain extent in the presence of the 
mixing. 

Now, we are in a position to calculate the critical temperatures Tc for the crossover transitions numerically. They 
are determined by the chiral and Polyakov susceptibilities as in Ref. 7J. Being almost equivalently, they can be also 
obtained from the maximum values of d^/dT and da^ jdT for the deconfinement and chiral phases, respectively as 



in Refs. [5 HTU1H5] . In the present talk, we employ the later method. In Table III, we list them for $ and tr^, assigned 
as Tf and T^, respectively, for each type of M. Note that we do not show the numerical results for for the cases 
with Mo,T and Mo,o, since they are UV divergent, proportional to / k^dk, unless a cutoff is introduced by hand. 

As discussed previously, from the table, one can see clearly that is lowered by inclusion of the (fc, T) dependence 
in M. It turns out that the shift of is about 10%, (240 227) MeV for Mq.q t- Interestingly, if we take into 

account fuU (fc, T) dependence, Tf and get closer to each other as shown in Table |lll| (T*, T'^') = (227, 216) MeV. 



About 5% discrepancy (^ 10 MeV) is observed between them, showing a tendency < T^- By turning off the T 
dependence in M {Al^fi), the discrepancy between and becomes larger up to about 15%, (T'^,T°' ) — (225, 265) 
MeV. For interpreting this behavior, we take a look on the T dependence of for the cases with 7- and M^fi, and 
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FIG. 4: $ (increasing curves) and (j^ jo^ (decreasing ones) as functions of T for To = 270 MeV, employing Affe^y (solid) and 
Mfc,o (long-dashed). 

draw the numerical results in Figure |4] For comparison, we also draw $ for the two cases. Note that show obvious 
difference for the cases with M^.t and Af^.o- This is the reason why the discrepancy between and is so large 
for Affc.Q. At the same time, this strong dependence on T for the chiral phase transition interprets the larger shift of 
T^, (265 — > 216) MeV. From this observation, we can conclude that the T dependence in M plays an important role 
in exploring chiral phase transition at finite T . 

It is worth mentioning that, from the LQCD analyses, it turned out that « 180 MeV for N f = 2 using the 
clover-improved Willson fcrmions Also, using the renormalization-group (RG) improved action, it was found that 
ss 171 MeV [5]. These values are significantly smaller than ours by (10 — 20)%. If this is the case, one may need 
more strong T dependence for M in the present approach, since the T dependence of is mainly governed by the 
behavior of M. In other words, the instanton effects must decrease much faster as T increases. There can be several 
possible scenarios to satisfy this condition: 

i) If we consider a correct instanton distribution function, rather than the simplified one, according to the large- A^^ 
limit, may be lowered. To test this, we have the ratio of M computed with the correct one: 



We, however, verify that this modification does not work for lowering T'^^ whereas the strength of Af is reduced 
approximately by half. 

ii) Additional T-dependent terms can be taken into account. For instance, we assumed that the Lgrangian multiplier 
A is independent on T in deriving M^^t- If it has the T dependence, we can modify Affc^T as follows: 



Quantitatively, the T-dependent A can not be determined self-consistently within the framework. From a very rough 
estimation, based on the assumption that no, which is the instanton packing fraction at T = 0, converges with a brute 
cutoff A « 1/po, we can obtain a relation Aq (x I/pq |33| . Using this assumption, the expression for M^^t is modified 
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TABLE III: Tc computed from $ and a^/ao for To = 270 MeV. 
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FIG. 5: $ as a function of T for To = 270. The solid and dashed Unes indicate the cases with and without the dynamical 
quark and $ mixing, respectively. The notations □ and A indicate the full and quenched lattice data, respectively, taken from 
Refs. gB] and [12]. 
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It turns out that this rough assumption makes things worse: is shifted to a larger value as expected in Eq. (171 



iii) We may not ignore the Matsubara frequency Wn in the denominator. In the presence of the mixing, $ may 
provide effects on the instanton distribution function, becoming an exponentially decreasing function, not a gaussian 
one. We, however, do not perform quantitative calculations for these possibilities here and leave them for future 
works. 

iv) Additional T dependence can be added to the instanton distribution function, according to the fermion overlap 
matrix, if one considers full QCD in computing the distribution function 43J, being different from the present talk 
based on the variational method in pure-glue QCD (the Harrington-Shcpard caloron) j31j . Moreover, the instanton 
clustering, which was suggested as a main contribution for the chiral phase transition [43l |44j and not taken into 
account here, may be responsible for lowering T^. 

v) We note that Tq for the pure-glue potential can be chosen as a smaller value than 270 MeV, which has been 
used throughout in the present talk. As shown in Refs. [5Ul0j. by taking Tq = 190 MeV, the computed values for 

became compatible with those from the LQCD simulations. Although we have not presented detailed results for 
lower To, we could obtain = 194 MeV using Mk,T at Tq = 200 MeV, showing about 10% decreasing. 

As for estimated in the LQCD, using the clover-improved Wilson fcrmions similarly, it was determined about 
210 MeV j6j, which is rather compatible with ours. It was observed in the usual local pNJL model [8nT0j that Tc 
taken from the two different order parameters are almost consistent: fa fa 220 MeV for Tq = 270 MeV. It 
turned out that = 215 MeV in Ref. with pNJL. 

Finally, we compare our results for $ with the LQCD data from Refs. [IB] and [13 in Figure [5] in which a full and 
quenched calculations were done for Nf = 2. In their works, it was observed that — 202 MeV, which is about 10% 
lower than ours, 227 MeV for Nf — 2, and 270 MeV for Nj = 0. As shown in the figure, the LQCD data, indicated by 
□ (full) and A (quenched), are in a qualitative agreement with the present results, but not quantitative. Especially, 
our result for the mixed case deviates much from it for the region T > Tq. 



IV. SUMMARY AND OUTLOOK 



In the present talk, we have attempted to derive an effective thermodynamic potential QoS at finite T and zero 
quark-chemical potential (/ip — 0) in the chiral limit. We restricted ourselves to A'c = 3 and Nf — 2. Motivated by 
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the Polyakov-loop-augmented Nambu-Jona-Lasinio model, we wanted to incorporate two different order parameters, 
and $, which characterize the chiral and deconfinment phase transitions, respectively. 

In order to discuss the spontaneous breakdown of chiral symmetry at finite T, we employed the singular-gauge 
instanton solution, and the fermionic Matsubara formula to express the effective chiral action as a function of T. We 
employed the instanton-distribution function, derived from the Harrington-Shepard caloron, to obtain the instanton 
density and average instanton size as functions of T. It turned out that these two quantities decreased but finite, 
indicating that the instanton effect survives even beyond [3T. The (fe, T)-dependent M, M^.t was derived in the 
large- iVc limit. We found that the k dependence of M becomes weaker as T increases. At the same time, the absolute 
value of Mfc T was also reduced with respect to T. To include the Polyakov loop as an order parameter for the 'L{Nc) 
symmetry, we took into account imaginary quark-chemical potential /x/ = ^4, which was indentified as the traced 
Polyakov loop <&. Combining all these ingredients, we could construct iloff with an additional pure-glue SU(iVt,) gauge 
potential. By minimizing rioS with respect to external fields such as a and <&, we could compute and as functions 
of T numerically. From the various numerical results we have found the followings: 

i) In the presence of the mixing of the dynamical quarks and we observed that $ is very sensitive to the mixing, 
showing the crosssover and first-order transitions with and without it, respectively. In contrast, is insensitive to 
it, indicating the crossover phase transition. 

ii) M was expressed as a decreasing function of T as well as k. Due to this, was lowered by about (5 — 10)%, 
in comparison to that with constant mass A/q.o ~ 350 MeV. From this observation, we explain this lowering Tf 
by that the nontrivial QCD vacuum contributions and partial chiral restoration play a significant role even in the 
deconfinment phase transition. 

iii) If the (fc,r) dependence had been fully taken into account, we found that Tf = 227 MeV and 216 MeV. 
The discrepancy between them became about 10 MeV, which was rather larger than that computed in the pNJL 
model. We also note that the LQCD simulations presented smaller than ours, whereas Tf was compatible. 

iv) Finally, we observed that was depending much on T. Again, the partial chiral restoration turned out to be 
crucial to make proper results for the chiral phase transition. 

Consequently, from the present talk, we could learn that the partial chiral restoration and nontrivial QCD vac- 
uum effects must be taken into account appropriately to investigate the breaking patterns of the chiral and Z(iVc) 
symmetries at finite T . As a next step, we attempt to include the finite quark-chemical potential (^r), giving a full 
description for the QCD phase diagram on the /iR-T plane. In addition to it, the finite current-quark mass, m is 
also under consideration beyond the chiral limit. This is important, since the explicit breakdown of the flavor SU(3) 
symmetry modifies the QCD phase diagram to a good extent. However, including finite m into the present framework 
has a huddle: One needs to consider the meson-loop corrections, which make significant modification on the physical 
quantities such as the chiral susceptibility, in comparison to those in the chiral limit [36] . Related works has been 
published [48j . and are under progress and appear elsewhere. 
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